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I. Introduction 

Spacecraft which utilize electric propulsion (EP) systems are capable of delivering a 
greater payload fraction compared to spacecraft using conventional chemical propulsion 
systems. Several researchers have investigated numerous applications of low-thrust EP 
including a manned Mars mission [1], scientific missions to the outer planets [2], and lunar 
missions [3]-[5] . In contrast, the study of optimal combined high and low-thrust spacecraft 
trajectories has been limited. 

In response to the release of NASA’s 1994 Announcement of Opportunity (AO) for 

Discovery class interplanetary exploration missions, a preliminary investigation of a lunar- 

comet rendezvous mission using a solar electric propulsion (SEP) spacecraft was performed. 

The Discovery mission (eventually named Diana) was envisioned to be a two-phase scientific 

exploration mission: the first phase involved exploration of the moon and second phase 

involved rendezvous with a comet. The initial phase began with a chemical propulsion 

translunar injection and chemical insertion into a lunar orbit, followed by a low-thrust 

SEP transfer to a circular, polar, low-lunar orbit (LLO). After scientific data was collected 

at the moon, the SEP spacecraft performed a spiral lunar escape maneuver to begin the 

interplanetary leg of the mission. After escape from the Earth-moon system, the SEP 
‘Final Report for NASA grant NAG3-1650 
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spacecraft maneuvered in interplanetary space and performed a rendezvous with a short- 
period comet. 

An initial study that demonstrated the feasibility of using EP for the lunar and comet 
orbit transfer was performed under the grant NAG3-1581 [6]. This final report is a contin- 
uation of the initial research efforts in support of the Discovery mission proposal that was 
submitted to NASA Headquarters in October 1994. Section II discusses the lunar orbit 
transfer phase of the Diana mission which involves both chemical and electric propulsion 
stages. Section III discusses the chemical lunar orbit insertion (LOI) burn optimization. 
Finally, section IV presents the conclusions of this research effort. 

II. Combined Chemical-Electric Propulsion Lunar Transfer 

The initial phase of the Diana mission involves a ballistic lunar orbit transfer, followed 
by a LOI chemical propulsion burn, and finally an EP orbit transfer to a polar, low-lunar 
100-km altitude orbit. The optimal lunar capture and circularization transfer using the 
solar electric propulsion (SEP) stage was outlined in the grant report NAG3-1581 [6]. In 
this section, the trajectory optimization study for the combined chemical-electric propulsion 
maneuver is presented. 

Trajectory Optimization 

The objective is to compute the minimum-fuel ballistic translunar trajectory from low- 
Earth orbit (LEO) to the optimal lunar orbit insertion (LOI) boundary conditions. The 
trajectory is shaped by two impulsive chemical burns at both ends. The complete optimal 
control problem is given below: 

For the free end-time problem, find the orientation and magnitude of the chemical burn- 
out velocity vector Vbo, the angular position of the spacecraft in LEO at translunar injection 
(TLI), and the magnitude of the chemical LOI AF which minimize 
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subject to the unpowered three-body equations of motion 
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with the initial conditions 

*(0) = <7(0o, 0o, v ho ) (8) 

and the terminal state constraints 

1 h p (tf) — h Vj \ / O'' 

SSf[X(tf),tj]= h a (tf) — h a{ = 0 (9) 

\i{t f )-9 0 deg/ [oj 

The motion of the coasting spacecraft is governed by the restricted three-body problem 
dynamics as indicated by Eqs. (2-7). The differential equations are formulated in a rotating 
Cartesian frame with the origin at the Earth-moon system center of mass and the positive 
x-axis pointing to the Earth along the Earth-moon line. The x-axis is fixed with the Earth- 
moon line, the y-axis is in the Earth-moon orbit plane, and the positive z-axis is along the 
angular momentum vector of the Earth-moon system. The Earth and moon are assumed 
to revolve in circular orbits about their common center of mass. The constant Earth and 
moon distances from the center of mass are denoted by x e and respectively. The 
constant angular rate of the Earth-moon system is u> and the gravitational parameters of 
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the Earth and moon are denoted by fi e and /x m , respectively. The position of the spacecraft 
in the rotating frame is denoted by ( x,y,z ) and the respective velocity components in the 
rotating frame are (u,v,w). The distances from the Earth and moon to the spacecraft are 
denoted by r e and r m , respectively. 

The initial conditions of the spacecraft at TLI are given by Eq. (8) and are a function 
of the position in LEO (do and i/> o) and the burn-out velocity Vbo as supplied by the upper 
stage of the Delta II launch vehicle. The initial velocity vq of the spacecraft with respect 
to an Earth-centered inertial frame is calculated by the below vector equation 

Vo - VLEO + Vbo (10) 

where vleo is the velocity of the spacecraft in a circular, 185-km altitude low-Earth orbit. 
By defining the longitude 0 o and heading tp 0 at TLI, the initial velocity veeo can be 
computed. The initial LEO is assumed to be inclined 7 deg with respect to the Earth- 
moon orbit plane which corresponds to launch conditions in late 2000. By specifying the 
launch energy C 3 and the pitch and yaw orientation of the excess burn-out velocity, the 
additional velocity Vb 0 can be computed. Finally, the position and velocity components in 
LEO must be transformed to the rotating frame centered at the Earth-moon system center 
of mass. 

The terminal state constraints as denoted by Eq. (9) require that the final lunar orbit 
after LOI be an elliptical polar orbit. Therefore, the first two terminal state constraints 
require perilune and apolune altitudes ( h p and h a ) match the desired values ( h Vj and h aj ) 
and the third constraint requires inclination i be 90 deg at the final time tj. 

Solution Approach 

Since our problem involves discrete control parameters, the optimal control problem is 
solved using a direct method. The optimal control problem is replaced with a nonlinear 
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programming problem (NLP) and the trajectory design variables are (7 3 , orientation of the 
spacecraft in LEO, orientation of vt , 0 , coast time tj, and magnitude of the LOI chemical 
AE. The NLP is numerically solved using sequential quadratic programming (SQP) which 
is a constrained parameter optimization method [7]. The SQP algorithm used here utilizes 
first-order finite differences to approximate the gradients and is due to Pouliot [8]. The SQP 
problem formulation involves seven optimization parameters and three equality constraints. 
The equations of motion are numerically integrated by using a standard fourth-order, fixed- 
step, Runge-Kutta integration scheme with 2000 steps. 

The performance index to be minimized is the negative mass after LOI. This is equiv- 
alent to maximizing m^oi- Spacecraft mass after LOI is computed by using the rocket 
equation: 


mLOi = mrue AV/c ( 11 ) 

where the AE is the velocity increment from the LOI burn and c = I sp g is the exhaust 
velocity of the chemical stage. Specific impulse I sp of the chemical stage is 310 s. The 
injected mass uitli is computed by using a simple linear fit of the launch performance of 
the Delta II: 

tutli = —27(73 T 1227 kg (12) 


Results 

Several optimal minimum-fuel lunar trajectories were readily obtained for a range of 
perilune and apolune altitudes by using the SQP optimization code. The optimal (7 3 was 
found to be —2 km 2 /s 2 for nearly every case which results in an injected mass of 1281 kg. 
Initially, the optimal circular LOI burn was obtained for a range of circular altitudes. That 
is, both h Pj and h af are equal and set at a wide range of altitudes from the lunar sphere of 
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influence (SOI) to the 100-km final low-lunar orbit (LLO). Since the optimal C 3 is nearly 
constant at —2 km 2 /s 2 , the impulsive AV required for the LOI burn has the most effect 
on rriLoi and the resulting optimal circular LOI Ak’s are presented in Fig. 1. The largest 
chemical propellant penalty is for a direct insertion into polar LLO (therefore bypassing 
the SEP transfer) as indicated by AV = 815 m/s in Fig. 1. The curve shows a minimum 
at a circular altitude of about 13,000 km where Ak = 580 m/s. 

Next, the optimal LOI burn for an elliptical lunar orbit was investigated. The desired 
apolune altitude h af was set at 50,000 km which is within the lunar SOI (which has an 
altitude of roughly 64,500 km). A range of perilune altitudes h Pf was utilized and the 
resulting optimal LOI AV is presented in Fig. 2. In this case, AV steadily decreases as 
perilune altitude is decreased and all cases showed better performance than the circular 
LOI maneuver. For safety purposes, the perilune altitude was chosen at 1000 km since 
trajectory simulations with a complete gravity model showed that a coasting trajectory 
with a perilune of 1000 km remained stable for several revolutions. The AV for this 
elliptical orbit is 223 m/s. Therefore, the 1000 x 50,000 km orbit is the starting point for 
the SEP circularization maneuver outlined in the previous final grant report [6]. 

The translunar trajectory to the 1000 x 50,000 km elliptical orbit is presented in Fig. 
3. Optimal transfer time is 5.6 days and the spacecraft “falls” into a vertical plane about 
the moon with a perilune altitude of 1000 km as shown by Figs. 4 and 5. Therefore, the 
optimal impulsive AV of 223 m/s is applied at perilune to produce the desired 1000 x 
50,000 km elliptical lunar orbit. 

Spacecraft Optimization with Chemical-EP Stages 

The optimal chemical translunar trajectory problem was extended to include spacecraft 
system optimization. In this problem, the chemical TLI and LOI burns were optimized 
along with the subsequent EP transfer to polar, circular LLO and the EP transfer to escape 
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Time ticks = 0.56 days 









Figure 5: Minimum-fuel translunar trajectory - above orbit plane 

conditions. The goal was to maximize the net mass of the spacecraft which is defined as the 
usable mass for payload plus basic spacecraft structural mass. Therefore, EP spacecraft 
system parameters such as I sp and input power P are optimized in order to determine total 
low-thrust propellant mass, EP tankage mass, and EP power and propulsion system mass. 
The quasi-circular transfers are calculated by using Edelbaum’s equations [9], The results 
were published in a refereed journal article [10] which has been attached to this final report. 

III. Optimal LOI Burn Maneuver 

Once the optimal translunar trajectory has been computed using impulsive AV com- 
putations, a realistic LOI burn optimization study is performed to determine the optimal 
steering profile for the chemical burn and the proper sizing of the insertion engines. The 
boundary conditions for the trajectory optimization problem are determined from the op- 
timal translunar trajectory from section II. 
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Trajectory Optimization 

The objective is to compute the minimum-fuel, finite-time chemical LOI burn which 
results in the desired 1000 x 50,000 km elliptical lunar orbit. A range of thruster levels 
was investigated to determine the best trade-off between thruster size and gravity loss. 

The complete optimal control problem is given below: 

For the free end-time problem, find the pitch and yaw thrust steering angles u(t) and 
v(t), and the burn duration 4 urn which minimize 


J = 

(13) 

subject to the three-body equations of motion 


dr 

di = Vr 

(14) 

dv r 

-7- = — + ar sin u cosu + vu T 
dt r 

(15) 

due v r v e 

—— — + ar cos u cos v + V Ug 

dt r 

(16) 

dfl sin# . . 

, = . . (o T sinu + Vt/fc) 

dt vo sin 1 

(17) 

di cos 6 . . 

= [a T sm v + V U k) 

(it Vq 

(18) 

d6 vo sin 6 cos i . . 

— = — — : — — far sin v + wUh) 

dt v vg sin 1 

(19) 

T 

where ar = r~ , 0 < t < t< 

tutli ~ 


with the initial conditions 


r(0) = 52,250.54 km 

(20) 

u r (0) = -0.921677 km/s 

(21) 

v e (0) = 0.108347 km/s 

(22) 
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ft(0) = 268.95 deg 


(23) 


z(0) = 89.18 deg (24) 

0(0) = 5.98 deg (25) 

and the terminal state constraints 

/ h p (tj)-mO km \ / 0 \ 

'h[x(ty),t/j = h a {tj) — 50,000 km = 0 (26) 

V *(</)- 90 deg ) V 0 / 

The states are radial position r , radial velocity u r , circumferential velocity longitude of 
the ascending node angle fi, inclination i, and in-plane longitude angle 6. The radius r is 
the distance from the center of the moon to the spacecraft and v r and vg are the inertial 
velocity components measured in the instantaneous orbit plane. The ascending node angle 
ft is measured counter-clockwise from the fixed +a: axis to the ascending node direction. 
The inertial +x axis is initially pointing from the moon’s center to the Earth at t = 0. The 
inclination i is with respect to the x — y or Earth-moon orbit plane. Longitude angle 6 is 
the in-plane angle measured from the ascending node to the spacecraft in the direction of 
motion. Therefore, 0 is the sum of argument of perilune u> and true anomaly i/. 

The gravity potential gradient VC/ for the combined Earth and moon gravity field is 


^ e ^ e—m 
/te ~ 


where the gravitational parameters of the Earth and Moon are represented by // e and // m , 
respectively, f e is the radius vector from the Earth to the spacecraft, r e _ m is the radius 
vector from the Earth to the moon, and D is the constant separation distance between the 
Earth and moon. The components of V(7 are 


W r = -^-^ + 


Vt/, = -^ + 

r e 


e—m r 

D 3 

(28) 

7’e—mg 

D 3 

(29) 
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V(4 = 


e/j l^e^'e—mh 

r 3 + D 3 ~ 


(30) 


where the subscripts r and d correspond to components along the radial and circumfer- 
ential in-plane directions and the subscript h corresponds to the direction normal to the 
instantaneous orbit plane. 

The in-plane pitch thrust steering angle u is measured positive above the local horizon 
to the projection of the thrust vector onto the orbit plane. The out-of-plane yaw thrust 
steering angle v is measured positive above the orbit plane to the thrust vector and is 
between ±90 degrees. The high-thrust acceleration of the spacecraft, aj, is computed by 
dividing the constant thrust magnitude, T, by the current spacecraft mass. Thrust levels 
ranging from 5 to 100 lb/ were investigated. The mass of the spacecraft is denoted by m, 
and propellant mass flow rate m is considered positive out of the vehicle. 

The initial conditions (20-25) represent a spacecraft state from the optimal translunar 
trajectory 13.7 hrs from perilune. At this given state, the spacecraft has crossed the SOI and 
has a moon-relative energy of 0.3368 km 2 /s 2 and eccentricity of 1.378. The three terminal 
state constraints (26) define the desired 1000 x 50,000 km polar elliptical orbit. The goal 
is to find the thrust steering angles u(t) and v(t), and the duration of the finite-burn arc 
such that the final spacecraft mass m(t/) is maximized (or, equivalently, such that fuel 
is minimized) and the spacecraft terminates in the prescribed polar elliptical orbit. The 
optimal control problem is solved using SQP with each continuous-thrust steering angle 
replaced by a cubic-spline fit through six SQP control parameters. 

Results 

Optimal minimum-fuel LOI burn maneuvers were readily obtained for constant thrust 
levels ranging from 5 to 100 lb/. The equivalent AV as computed by the rocket equation 
for the resulting finite burn is presented in Fig. 6. The performance greatly improves and 
asymptotically approaches the impulsive AV of 223 m/s as thrust level is increased. The 
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gravity loss can be computed as the difference between the finite AV and the impulsive 
AV. The “knee” of the curve is at a thrust range of about 20-35 lb/ where the gravity 
loss ranges from about 18 to 8 m/s. The optimal finite-burn duration is presented in Fig. 
7 and a similar profile is observed with the burn duration ranging from 54 to 25 min for 
the “knee” of the curve. Using these plots, it was determined that a thrust level of 24 lb/ 
would provide sufficient performance. This thrust level could be easily supplied by three 
8-lb/ thrusters that are currently available from TRW. 

The resulting optimal pitch steering history for the 24-lb/ thruster is presented in Fig. 
8. The pitch steering angle plotted here is measured from the velocity vector and shows a 
linear relation symmetric about the 180 deg steering angle. Therefore, the optimal pitch 
steering profile initially points opposite the velocity vector but with a slight positive radial 
velocity component. The thrust vector is then turned at a constant rate in the orbit plane 
and is directly opposite the velocity vector (180 deg) at the mid-point of the burn. The 
optimal yaw steering is effectively zero since the maneuver is essentially planar. The burn 
lasts for about 45 min. The optimal LOI perilune burn is presented in Fig. 9 as viewed 
from the moon-centered, inertial y-z vertical plane. 

IV. Summary and Conclusions 

A study of a translunar injection trajectory using a combined chemical-electric propul- 
sion spacecraft has been performed. The study was in support of the Diana mission proposal 
for NASA’s 1994 Announcement of Opportunity (AO) for Discovery class exploration mis- 
sions. The trajectory optimization was performed using sequential quadratic programming 
(SQP), which is a direct method. 

The optimal ballistic translunar trajectory was obtained using SQP and impulsive burns. 
A three-body dynamic model was used for the governing equations of motion. A complete 
vehicle sizing study for the combined chemical-electric propulsion spacecraft was performed 


13 



Burn Duration , min. 



Figure 6: Optimal AV vs chemical thrust level 



Figure 7: Optimal burn duration vs chemical thrust level 
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and the results were published in a refereed journal. 

A detailed study optimizing the lunar orbit insertion burn was performed. A finite- 
duration perilune burn was used and a range of thrust levels was investigated in the context 
of the complete three-body equations of motion. A good trade-off between thruster size and 
gravity loss was identified and a 24-lb/ thruster was selected. The optimal thrust steering 
history was obtained and the result is a simple linear pitch profile. 
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